Rate Equations and Scaling in Pulsed Laser Deposition 
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We study a simplified model for pulsed laser deposition [Phys. Rev. Lett. 87, 135701 (2001)] 
by rate equations. We consider a set of equations, where islands are assumed to be point-like, as 
well as an improved one that takes the size of the islands into account. The first set of equations is 
solved exactly but its predictive power is restricted to a few pulses. The improved set of equations 
is integrated numerically, is in excellent agreement with simulations, and fully accounts for the 
crossover from continuous to pulsed deposition. Moreover, we analyze the scaling of the nucleation 
density and show numerical results indicating that a previously observed logarithmic scaling does 
not apply. 



PACS numbers: 68.43.Jk, 



1.55.A-, 05.45.Df 



I. INTRODUCTION 

Among the techniques for the growth of thin films and 
multilayers Molecular Beam Epitaxy (MBE) and Pulsed 
Laser Deposition (PLD) [1] play an important role. Com- 
pared to MBE, PLD has several advantages. For exam- 
ple, pulsed laser deposition is a widely used technique 
which allows one to control the stoichiometry of multi- 
layers more efficiently and leads in some situations to a 
better layer-by-layer growth [2| . The essential difference 
between these two techniques lies in the way how parti- 
cles are deposited and the kinetic energy of the arriving 
particles. While in MBE the flux of particles is contin- 
uous, in PLD intense and short pulses are deposited on 
the target, making the flux of incoming particles approxi- 
mately discontinuous. Also, in PLD the kinetic energy of 
the deposited particles is typically higher than in MBE. 

There is already a great amount of theoretical work 
on MBE [aL while there are less theoretical studies on 
PLD 0, [ill 0, & S G3- In the present paper we ap- 
ply to PLD a theoretical method that proved to be very 
useful in MBE, namely, the rate equations approach [llj | . 
By integrating a set of rate equations for the island den- 
sity one can make predictions about the experimentally 
relevant quantities characterizing the system. We use 
two sets of rate equations: In the first and more simple 
one, we consider the islands as point-like objects, while 
in the second approach, introduced for MBE in Ref. [l^] , 
we improve the equations by considering the competition 
of islands of different sizes for the arriving and diffusing 
monomers. The first approach has the advantage that 
the rate equations can be solved exactly. The improved 
set of rate equations can only be integrated numerically 
but it shows excellent agreement with simulations. 

In this work we study PLD on the basis of a model 
investigated in Ref. where monomers are deposited 
on the surface in infinitely short pulses. These monomers 
diffuse on the surface between subsequent pulses until 
they nucleate to form the seeds of immobile islands. More 
precisely, the model is defined as follows: 



- In each pulse I ■ L? monomers, where I is the 
pulse intensity, are randomly deposited on a two- 
dimensional L x L lattice. 

- During the time interval between two pulses of 
length t = I/F, where F is the time-averaged in- 
flux of particles, the monomers diffuse with rate D. 

- When a monomer encounters another monomer or 
the border of an agglomerate of particles, it sticks 
irreversibly. 

The model does not take into account that the arriving 
particles in PLD have a high kinetic energy, leading to 
transient effects such as an enhanced mobility. In fact, 
the transient mobility was shown to have a great influ- 
ence on the growth kinetics [j| H [l(|. However, Vasco 
et al. [ioj . analyzing a set of rate equations, found out 
that only the way particles are deposited already leads 
to differences in the film roughness. More specifically, 
they found that continuous deposition leads generally 
to rougher films and an additional transient mobility in 
the case of continuous (discontinuous) deposition makes 
the film rougher (smoother). Other simplifications as- 
sumed in the model are the irreversibility of aggrega- 
tion and the immobility of nucleated islands. For in- 
stance, the assumption of irreversible aggregation hinders 
to study Ostwald ripening, which was recently proposed 
to be the key mechanism ruling the PLD growth kinet- 
ics in momoner-depeted regimes (Toj . Here we want to 
study the difference in the island morphology for contin- 
uous and discontinuous deposition and the correspond- 
ing crossover in the limit of low energies, assuming that 
transient mobility effects lead to corrections that can be 
neglected. Moreover, we restrict ourselves to the sub- 
monolayer regime, where the coverage 0, which is given 
by the pulse intensity multiplied by the number of pulses, 
is less than one. It turns out that the dynamical processes 
at the bottom layer determine the _growth process at the 
following layers to a large extent [13] ■ 

A nucleation event occurs when a diffusing monomer 
encounters another monomer. A quantity of interest in 
the study of PLD and MBE is the nucleation density n, 
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which is defined as the number of nucleation events di- 
vided by the number of lattice sites. In Ref. [f| it was 
suggested that for PLD this quantity follows an unusual 
logarithmic scaling. Later this type of scaling behav- 
ior was explained in the framework of scaling laws with 
continuously varying exponents [14j . Opposed to these 
early conjectures, we present numerical evidence that 
such logarithmic scaling laws do not hold asymptotically, 
although they may be used as a good approximation. 

The paper is organized in the following way. In the 
next section we present the properties of the model and 
introduce the observables that we are going to calcu- 
late. The third section is dedicated to the rate equations, 
where islands are approximated as point-like objects. In 
section 4 we introduce and analyze the improved set of 
rate equations. Section 5 is concerned with a critical 
analysis of scaling laws for PLD. Finally, our conclusions 
are summarized in section 6. 



II. MODEL PROPERTIES 

Depending on the pulse intensity / and on the ratio 
R = D / F the model for PLD defined above displays 
different features. For very small pulse intensities and a 
finite value of R, it exhibits essentially the same behavior 
as MBE. More specifically, for intensities much smaller 
than a typical value I c , which depends on R, the model is 
in the MBE regime. Then, increasing the pulse intensity 
to a value much larger than I c , the model crosses over 
to a different regime, the so-called PLD regime. The two 
regimes differ in so far as they are characterized by differ- 
ent surface morphologies, that is, the number and size of 
islands are distributed differently. However, as edge dif- 
fusion is not included in the model, the islands remind of 
fractal objects in both cases and their fractal dimension 
df is approximately equal to the one of diffusion limited 
aggregation (DLA) [TEGl]. [H 

The two regimes are also different with respect to the 
scaling of the average distance between islands Id ■ In the 
MBE limit / <C I c the average distance between islands 
grows as Q 



Relations (QJ and ([3]) can be combined in the following 
scaling form Q 



where 



In ~ i? 7 , 



1 



d,+4' 

while in the PLD limit I ^> I c one has [|[ 
where 

i/ = 7 /(l -27). 
Hence, the crossover takes place if 



(1) 
(2) 

(3) 
(4) 
(5) 



l D ~ R 1 h{I/R- l/v )i 



(6) 



where h(x) is a scaling function which is constant for 
1 < 1 and h(x) ~ x~ v for x ^> 1. The same type of 
scaling behavior is also expected to hold when consider- 
ing other typical length scales as, for example, the mean 
distance that a monomer travels before being captured 
by an island or another monomer, the square root of the 
mean island size, and the inverse of the square root of the 
nucleation density. In numerical simulations the scaling 
relations |T]) , ^ and © were observed to hold for small 
coverages before coalescence of islands sets in (see below) . 

In the PLD limit the model is independent of R, the 
time interval between two pulses is big enough to make 
the density of monomers equal to zero. Hence, taking the 
PLD limit means that diffusion of monomers occurs until 
the density of them is zero, and in principle such a limit 
can be taken for any value of /. The smallest value of / 
with which we perform calculations is 0.0001. This value 
is much lower than typical experimental pulse intensities, 
which are of order 0.1 [H,[ll. We use low values for the 
pulse intensity because we are interested in the scaling 
behavior of the model in the limit / — * 0. 

The probability that an occupied site belongs to a clus- 
ter of size s is given by 



Ps 



(7) 



where N s is the number-density of islands of size s and 

00 

9 = Ft = sN s (8) 

s=l 

is the coverage. Therefore, the mean island size is given 
by 



(9) 



The nucleation density n is defined as the number 
of nucleation events (when a monomer meets another 
monomer) in the bottom layer per unit area. Initially, 
when freshly nucleated islands do not yet touch each 
other, it is equal to the island density, given by 



n = J^n s . 



(10) 



However, as the coverage increases, more and more is- 
lands coalesce, forming bigger islands. Consequently 
their number decreases so that the island density be- 
comes smaller than the nucleation density, as is shown 
in FigfTJ As can be seen, coalescence - indicated by a 
difference between n and N - sets in already before the 
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up by a theoretical framework which involves continuous 
varying exponents 14|, where the scaling exponents are 
continuous functions of the control parameters. One of 
the purposes of this paper is to show that this logarith- 
mic scaling does probably not apply to PLD in a strict 
sense. 

Throughout this paper all numerical simulations were 
performed in 2+1 dimensions on a 400 x 400 lattice with 
periodic boundary conditions. For average quantities 
such as the island and nucleation densities or the mean is- 
land size, the number of independent runs was 100 while, 
when calculating the probability distribution ((7|) , we used 
10000 runs. 



FIG. 1: (Color online) Nucleation density n(f) and island 
density N(t) obtained from simulations in the PLD limit for 
I = 0.001. 



island density reaches its maximum. At this maximum 
the island density begins to decrease because coalescence 
events become more frequent than nucleation events. 

As for the scaling of the nucleation density in the PLD 
regime, the following relations were observed in Ref. ||: 
After the first pulse, where = 1, the nucleation density 
scales as 



n(/,/)~7. 

Upon completion of the bottom layer 
ation density was found to scale as 

n(I,l)~I 2 ». 



(11) 

1) the nucle- 



(12) 



These two scaling laws determine the terminal points of 
the curves n(I, 9) for different intensities /. However, 
these curves were found to have a different shape in a 
double-logarithmic representation, making it impossible 
to perform an ordinary data collapse. This suggested 
that the model does not exhibit ordinary power-law scal- 
ing. 

As a possible way out, a logarithmic scaling form was 
proposed Starting point was the observation that 
by squeezing and stretching the curves in such a way 
that the terminal points collapse, the entire curves for 
different intensities seem to collapse as well. This led to 
the conjectured scaling form 



lnM(I,9) = QnI)g(^ 



where 



M(I,6) = 



n(I,l) 



(13) 



(14) 



Here g(x) is a universal function, and this scaling is sup- 
posed to be valid in the full range < 6 < 1. In a sub- 
sequent work, this type of scaling behavior was backed 



III. RATE EQUATIONS 

In this and the following section we study the model 
introduced above by mean-field approximations on dif- 
ferent levels. Before starting let us emphasize that these 
approaches are based on the same assumptions (e.g. ir- 
reversible aggregation, immobile aggregates) as the full 
model defined in the introduction, extended by additional 
mean-field type approximations such as the homogeneity 
of densities and the neglect of correlations. 

We start with a simple set of rate equations for the 
island density ./V and the monomer density N% (see 
Ref. @,E[I|) 



dNi 



F(t) — 2DNi — DNNi, 



dN 
~dt 



DNi 



(15) 



(16) 



The first term on the right hand side of Eq. (fT5"|) de- 
scribes the flux of particles while the second and third 
terms are related to nucleation of two monomers and at- 
tachment of monomers at existing islands, respectively. 
In Eq. (|16p the term on the right hand side describes 
the increase of the island density by nucleation of two 
diffusing monomers. For MBE the flux is constant, while 
for PLD the flux is modeled by a discontinuous sequence 
of spikes 



F(t)=lJ2S(t-rl), 



(17) 



2=0 



where r is the time interval between two pulses. 

In this set of rate equations coalescence is not taken 
into account, hence within this theory the nucleation den- 
sity and the island density are identical. Moreover, the 
approach does not take into account how islands of differ- 
ent sizes compete for the diffusing and arriving monomers 
since islands are treated as point-like objects. Therefore, 
this approximation is expected to be valid (when com- 
pared with simulations) just for the first few pulses. 
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The MBE limit corresponds to taking r — » and / — ► 
with J/t = i 7 ", since this renders a constant flux F. 
In the MBE regime it was found that initially, when 
Ni ^> N, the island density grows as iV ~ t 3 . How- 
ever, when N ^> Ni, the growth of the island density 
becomes slower, crossing over to the power law N ~ i 1 / 3 
(see Ref. @). 

To solve the equations in the PLD limit, let us first 
consider the evolution between two pulses / and I + 1 : 



™L = -2N?-NN 1 , 
dt 1 



dN 
~dt 



= Nf, 



(18) 



(19) 



where, without loss of generality, we set D = 1. In or- 
der to make the above equations linear one introduces a 
modified time variable [HI, [22j 



T : 



N^dt'. 



(20) 



Because of dT = Ni(t)dt, Eqs. (JISJ) and ([11]) turn into 
dNx 



dT 



dN_ 
~dT 



-2N X - N, 



(21) 



(22) 



With the initial conditions Ni(0) = I and N(0) = N®, 
where is the island density after I pulses, one is led 
to the solutions 



Ni(T) = [I-(I+N^)T] exp(-T), 



(23) 



N{T) = [N (i) + + N®)T] exp(-T). (24) 

In the PLD limit the pulses are so strongly separated 
in time that all diffusing adatoms have nucleated or at- 
tached to existing islands before the next pulse arrives. 
The diffusion process ends at the final (modified) time 
Tf when 



Ni(T f )=0 



(25) 



and the corresponding island density N(Tf) = N^ l+1 ^ 
stays constant until the next pulse arrives. Because of 
Eq. (f2"3"| the modified final time Tf is given by Tf — 
I/(I+N®), hence with Eq. (24|) we obtain the recurrence 
relation 



N (W) = (j + jv(0) 



exp 



-/ 



(26) 



with the initial condition A^ ) = 0. Rewriting the island 
density as 



= Ail 



(27) 
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FIG. 2: (Color online) Island density in the PLD limit as a 
function of I for the first three pulses Z = 1,2,3 obtained from 
simulations (points) and from relations (|27[1 . (|28p . and (I29p 
(lines). 



with 



this recurrence relation turns into 



.4 



l+i 



(1 + Ai) exp 



l + A, 



(28) 



(29) 



which is independent of /. This means that within this 
mean field approximation the island density is propor- 
tional to the pulse intensity. As shown in Fig. [21 where 
relations (|2"T)) . |28p and (f2"9"]l are compared with simula- 
tions, this is indeed the case for the first few pulses. For 
the predicted slope, however, we see a better agreement 
for the first pulse and also for smaller values of the pulse 
intensity. 

This approximation reproduces another important fea- 
ture of the full model, namely, after the first pulse, for a 
fixed value of J, the nucleation density in the PLD limit 
is considerably larger than the nucleation density in the 
MBE limit R < /V^t-i) ( see Fig. H where the rate 
equations and simulation results are compared in both 
limits). However, the nucleation density in the MBE 
limit grows faster and eventually exceeds the nucleation 
density in the PLD limit such that the two curves in 
Fig- LH cross each other. 

The rate equations used here assume islands to be 
point-like objects. Obviously, this approach works only 
for very small values of the coverage (much smaller than 
the value at which coalescence starts) and thus it is in- 
adequate to fully account for the crossover from PLD 
to MBE. In particular, it does not allow one to confirm 
the scaling relation ([5]). In the next section we consider 
improved rate equations that overcome these problems. 
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FIG. 3: (Color online) Comparison of the nucleation density 
obtained by simulations (points) and the rate equations (|15[) 
and dT6l) (lines) with I = 0.001 and, in the MBE limit, R = 
10 4 . Note that for comparing the two cases for given 7, R has 
to be chosen such that R < jVC*?-*), 



IV. IMPROVED RATE EQUATIONS 

In order to improve the rate equations for the model 
defined in the introduction let us now consider the dy- 
namics of islands with different sizes separately and study 
how they compete for the diffusing monomers and how 
other adatoms are deposited on top of these islands. 

A set of rate equations for MBE that takes these fea- 
tures into account is [121 



^ = (D/F)Nx(a a - l N a - l -a a N a ) 
+fc s _iiV s _i - k s N s , 



(30) 



7AT OO 

OO 

-feiJVi -£fc s A s , (31) 

s=l 

where the constants a s are effective rates describing how 
an island of size s competes for the diffusing monomers 
while the constants k s are related to the capture of the 
incident monomers. These quantities are given by [l2| 



K = s 2/df 



and 



Ki(Rs/0 
K (R s /0 



(32) 



(33) 



where K and Ki are modified Bessel functions, R s is the 
island radius, and £ is the average distance a monomer 



travels before being captured by an island or another 
monomer, given by 



(34) 



s=2 



These improved relations assume that the islands are cir- 
cular. However, in reality and simulations the islands 
resemble fractal objects. The fractal morphology of the 
islands is taken into account in an effective way by as- 
suming that the island radius grows with its size by [e| 



R s 



a s 



l/df 



(35) 



where a is a fit parameter, to be determined by the com- 
parison of the results arising from numerical integration 
of the rate equations and simulations, and df — 1.72 is 
the fractal dimension of DLA-like clusters [23| . 

In order to generalize the above rate equations to 
pulsed deposition we use, as in the former case, a discon- 
tinuous flux of incoming monomers, obtaining the differ- 
ential equations 



dZVi 



£(l 16(0 - II) 

I ^ 8=1 ' 

oo 

-2Da x N\ - DNi £ <r s N s , (36) 



dN a 
d9 



Y^[K-iN s -i - k s N s ]IS(6 - II) 
i 

+DN 1 [a s _ 1 N s _ 1 ~ ( j s N s ] 1 (37) 



where the ^-functions couple to the terms proportional to 
k s since those terms are related to the capture of arriv- 
ing monomers, an event that takes place during deposi- 
tion. As in the previous case, the improved equations 
predict that the nucleation density n and island den- 
sity ./V are identical since coalescence is not taken into 
account. Since coalescence influences directly only the 
island density, we find it more appropriate to compare 
the results coming from numerical integration of these 
equations with the nucleation density obtained by simu- 
lations. 



The rate equations ([36]) and (|37|) faithfully reproduce 
the predicted crossover from MBE to PLD. To demon- 
strate this we show, in Fig. [51 the corresponding data 
collapse according to the scaling form © with v — 0.23 
and 7 = 0.15, where we used a = 1. We note that the 
value of the fit parameter a could be different for MBE 
and PLD. However, when considering the crossover from 
MBE to PLD one has to use a single value for which 
a = 1 turned out to be a good choice. This approxima- 
tion is justified since the value of this parameter does not 
change the type of scaling behavior. Nevertheless, this 
ambiguity regarding the value of a may lead to small 
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FIG. 4: (Color online) Data collapse to verify equation © FI G- 6: (Color online) The probability distribution function 
obtained from numerical integration of (J36J) and @2J with © for I = 0.0001 in the PLD limit obtained from simulations 
a — i a t tjjg coverage 9 = 0.1. (S) an d from a numerical integration of the rate equations (R) 

£M- and dHTJ) . 
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FIG. 5: (Color online) Comparison between simulations (S) 
and numerical integration of the rate equations (|36[1 and (|37[) 
(R) of the nucleation density (upper panel) and mean island 
size (lower panel) in the PLD limit for the pulse intensities 
/ = 0.0001 (a), / = 0.0002 (b), I = 0.0005 (c), and / = 0.001 
(d). 



numerical deviations in the estimates of the scaling ex- 
ponents v and 7. 

In Figs. [5]and[6]we compare the numerical integration 
of the rate equations with numerical simulations of the 
full model. Fig. [5] shows the nucleation density and the 
mean island size as functions of the coverage for different 
values of the pulse intensity. As can be seen, simulations 
and numerical integration of the rate equations agree al- 
most perfectly. As is shown in Fig. [SJ the agreement is 
less good for the probability distribution ([6]). Such dis- 
crepancies are not surprising because the rate equations 
approach is a mean field theory, neglecting correlations in 
the system. Similar discrepancies were already observed 



in the MBE limit [l2j]. The values of the fit parameter 
that we used are a « 1.7 in the PLD limit and a ~ 0.3 
in the MBE limit. 

This improved theory has considerable advantages 
when compared to the previous one. As shown in Fig. |4l 
it fully accounts for the crossover from MBE to PLD, 
and the observables, except for coverages where coales- 
cence starts to play an important role, are in excellent 
agreement with simulations. This excellent agreement in- 
dicates that the improved mean field equations capture 
the most important features of the model while other 
possible sources of deviations such as correlations and 
fluctuations are probably less important. 



V. SCALING PROPERTIES OF PULSED 
DEPOSITION 

As mentioned in the introduction, it was proposed that 
PLD could be described by an unconventional logarith- 
mic scaling form [5j. In this section we study this con- 
jecture from a critical perspective. 

Let us first recall the scaling properties of MBE. After 
an initial transient time and before the onset of coales- 
cence, the probability distribution ([7]) is known to obey 
the scaling relation [12| 

Ps(0) = (s)- 1 f^, (38) 

where f{x) is a universal function independent of the 
control parameters R and 9. This means that, in the 
MBE regime, the model has only one characteristic size 
and therefore the morphology of the islands stay the same 
if the control parameters are varied. 
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FIG. 7: (Color online) The probability distribution p a , ob- 
tained from simulations, multiplied by mean island size (s) as 
a function of s/{s) for I = 0.0001 and different values of 9 
(left); for 9 = 0.1 and different values of I (right). 



Turning to the PLD limit, in Fig. [7] we observe that a 
convincing data collapse is obtained for different cover- 
ages if the pulse intensity is kept fixed, while the collapse 
is not satisfactory for different values of the pulse inten- 
sity and a fixed coverage. Therefore, in the PLD limit, 
the function f(x) is not universal, instead it explicitly 
depends on one of the control parameters, namely, the 
pulse intensity I. However, for reasons that are explained 
below, the island morphology seems to be the same for 
different values of the pulse intensity. 

To explain the failure of conventional scaling, the log- 
arithmic scaling form (|T3|) was suggested on a purely 
heuristical basis This scaling form had been ap- 

plied previously to various other problems: experiments 
in turbulence [2J] , self-organized critical sandpile models 
[IS HH, and, as in the present case, DLA-related growth 
processes [27| . Fig. [5] shows a data collapse of the nor- 
malized nucleation density according to this scaling form, 
which at first glance seems to be convincing. Initially it 
was speculated that the unusual type of scaling behavior, 
which can be explained in terms of continuously vary- 
ing critical exponents |l4|, may be related to the fractal 
structure of the islands, which become more and more 
compact as the first monolayer is filled up. However, 
later it was shown [7J] that the same logarithmic scaling 
form can be used in a 1+1-dimensional model for pulsed 
deposition, where the (one-dimensional) islands are al- 
ways compact. 

A closer look at Fig. [8] reveals that the collapse is not 
perfect, rather there are deviations with a clear system- 
atic tendency. It seems that the curves become more 
straight as the intensity is reduced, questioning the con- 
cept of logarithmic scaling. 

Although we are still unable to disprove or confirm the 
concept of logarithmic scaling applied to PLD, it is in our 
opinion useful to demonstrate that an alternative scaling 
concept yields at least as good if not even better results. 
Starting point is the observation that average quantities 
like the mean island size and the nucleation density ex- 
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FIG. 8: (Color online) Data collapse of the normalized nu- 
cleation density, coming from simulations data, with the log- 
arithmic scaling form. The inset shows the interval 0.2 < 
In 0/ In I < 0.25 where the scaling function increases with the 
pulse intensity I. 



n(Q,I) 




FIG. 9: (Color online) Nucleation density, obtained from sim- 
ulations, as a function of I for different values of 9. The graph 
suggests that the exponent a increases with 9. 



hibit a clean power-law dependence on / if the coverage 
9 is kept fixed. However, the value of the exponent varies 
with the coverage, i.e. 



i(I,6) ~ I a{e) . 



(39) 



The plots of the nucleation density as a function of / 
for different values of 9 (cf. Fig. [9]) suggest that the 
exponent a increases monotonically with 9. Actually, 
a = 2v and therefore Eq. ((3]) should be corrected by 
introducing an exponent v which depends continuously 
on the coverage. Since the exponent v is related to the 
fractal dimensionality of the islands, this indicates that 
the fractal dimension of them effectively varies with 9. 
However, the data collapse of the probability distribution 
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for different values of the coverage (Fig. [7]) shows that 
this is clearly not the case. 

We have also observed that in the one dimensional 
model, where the islands have no fractal properties, the 
nucleation density follows a power law for a fixed value 
of the coverage. Moreover, we observed that the function 
f(x), defined in Eq. ([55)1 . does depend on the pulse inten- 
sity. This indicates that these two features of the model 
in the PLD limit are not related to the fractal dimension 
of the islands and hence we do not expect it to vary with 
the control parameters. 

In Fig. [8] we present a data collapse using relation (fl3|) 
for different values of the pulse intensity. Since M(7, 1) = 
1 and M{I, I) follows a power-law, the first (9 = 7) and 
last (9 = 1) points of the curves collapse. However, the 
situation in the middle of the curve (7 <C 9 -C 1) is 
unclear. More specifically, as is shown in the inset of 
Fig. |U for 0.2 < In 6/ In I < 0.25 the value of lnM(7,0) 
increases monotonically with 7, thus exhibiting a system- 
atic deviation rather than a statistical error. Also, this 
scaling is in disagreement with relation (I39|) . since from 
it the value of M(7, 9) should be constant for a fixed 
coverage. For these reasons we consider the numerical 
evidence supporting equation (|39[) as more reliable than 
the data collapse obtained using the logarithmic scaling 
form. 



VI. CONCLUSION 

In this paper we have studied several variants of rate 
equations for pulsed deposition. First a simple set of 
rate equations, where the islands are treated as point-like 
objects, was considered. The equations were solved ex- 
actly in the limit of strong and temporally separate pulses 
(the so-called PLD limit as opposed to the MBE limit of 
continuous deposition) , reproducing some features of the 
model for very small values of the coverage. Since this 
approximation does not take the dimension of the islands 
into account it was not possible to faithfully reproduce 
the crossover from MBE to PLD. This problem was over- 
come with a second improved set of rate equations that 
takes the dimension of the islands into account. As in the 
case of MBE [12| , we showed that these improved equa- 
tions lead to results which are in excellent agreement with 
simulations for pulsed deposition. 

Another point of the present work was to revisit a re- 



cently proposed logarithmic scaling for pulsed deposition 
from a critical perspective. In the corresponding data col- 
lapses we have observed small violations for intermediate 
coverages with a systematic drift, indicating that loga- 
rithmic scaling may be a good approximation for pulse 
intensities in computer simulations and experiments but 
probably it is not asymptotically valid in the limit 7 — > 0. 

As an alternative suggestion, we have proposed that 
the nucleation density, for a fixed value of the coverage, as 
a function of 7 follows a power law and that the exponent 
varies with the coverage. This suggestion is supported 
by numerical simulations and leads to numerical results 
which are more accurate than the data collapses obtained 
by using logarithmic scaling. This may be another hint 
that logarithmic scaling, as proposed in Ref. Q, has to 
be replaced by a different type of scaling theory. 

We have pointed out another new feature of the two- 
dimensional model in the PLD limit, namely, the proba- 
bility distribution p s that a site belongs to a cluster of size 
s does not scale according to relation (f3"5|) for different 
values of the pulse intensity. At first glance this non- 
universal behavior may lead to the conclusion that for 
different pulse intensities the fractal dimension of the is- 
lands is different. However, studying the one-dimensional 
model, where islands have no fractal properties, the same 
effect is observed, indicating that this is not the case. 

Let us finally comment on the one-dimensional model 
in more detail. This case was previously studied in 
Ref. 0|, where a reasonable data collapse for the 
nucleation density based on logarithmic scaling was 
presented. Contrarily, we observe a clean power law 
behavior of the nucleation density as a function of 7 for a 
fixed 9 and also a non-universal probability distribution 
for different values of 7. Hence, the one dimensional 
model seems to exhibit a behavior which is as rich as in 
the two-dimensional case and a systematic study of the 
Id case may be very useful for clarifying some aspects of 
PLD. 
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